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Observational results of compact objects are best understood using advec- 
tive accretion flows (Chakrabarti, 1996, 1997). We present here the results of 
numerical simulations of all possible types of such flows. 

Two parameter (specific energy £ and specific angular momentum A) space 
of solutions of inviscid advective flow is classified into 'SA' (shocks in accretion), 
'NSA' (no shock in accretion), T (inner sonic point only), 'O' (outer sonic point 
only) etc. (Fig. 1 of Chakrabarti, 1997 and references therein). Fig. la shows 
examples of solutions (Molteni, Ryu & Chakrabarti, 1996; Eggum, in prepara- 
tion) from 'SA', T and 'O' regions where we superpose analytical (solid) and 
numerical simulations (short dashed curve is with SPH code and medium dashed 
curve is with TVD code; very long dashed curve is with explicit /implicit code). 
The agreement is excellent. In presence of cooling effects, shocks from 'SA' os- 
cillate (Fig. lb) when the cooling timescale roughly agrees with postshock infall 
time scale (Molteni, Sponholz & Chakrabarti, 1996). The solid, long dashed 
and short dashed curves are drawn for T 1 / 2 (bremsstrahlung) , T 0A and T 0,75 



Fig. 1: (a) Numerical simulations and theoretical calculations of shock- free (O, I) and 
shocked (SA) flows, (b) Oscillation of shocks when cooling effects are included. 



cooling laws respectively. In the absence of steady shock solutions, shocks for 
parameters from 'NSA' oscillate (Fig. 2) even in the absence of viscosity (Ryu 
et al. 1997). The oscillation frequency and amplitude roughly agree with those 
of quasi-periodic oscillation of black hole candidates. When the flow starts from 
a cool Keplerian disk, it simply becomes sub-Keplerian before it enters through 
the horizon. Fig. 3a shows this behaviour where the ratio of \/\Keplerian is 
plotted. When the flow deviates from a hot Keplerian disk, it may develop a 
standing shock as well (Fig. 3b) (Molteni et al. 1996). 



Fig. 2: Oscillation of shock location 

when the flow has two sonic points 

but a steady shock condition is not satisfied. 



Fig. 3: (a) Ratio Xdisk/^Kep as obtained from a ID numerical simulation, (b) Another 
ID simulation which shows formation of sub-Keplerian flows and shocks from a hot 
Keplerian disk. 
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